Percolation and Magnetization 
in the Continuous Spin Ising Model 
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Abstract. In the strong coupling limit the partition function of SU{2) gauge 
theory can be reduced to that of the continuous spin Ising model with nearest 
neighbour pair-interactions. The random cluster representation of the continu- 
ous spin Ising model in two dimensions is derived through a Fortuin-Kasteleyn 
^ ! transformation, and the properties of the corresponding cluster distribution are 

Q I analyzed. It is shown that for this model, the magnetic transition is equivalent to 

O ■ the percolation transition of Fortuin-Kasteleyn clusters, using local bond weights. 

These results are also illustrated by means of numerical simulations. 

1 Introduction 



It has recently been proposed that the deconfinement transition in SU(2) gauge theory can be 
characterized as percolation of Polyakov loop clusters |l| . This idea is based on the fact that 
SU(2) gauge theory and the Ising model belong to the same universality class |0, and that 
the magnetization transition in the Ising model can be specified as percolation of clusters 
defined through local bond weights 0. It thus seems natural that the same holds for the 
corresponding Polyakov loop clusters in SU(2) gauge theory, provided suitable bond weights 
can be defined. In [|l|, it was shown that in a lattice formulation effectively corresponding 
to the strong coupling limit, SU(2) gauge theory indeed leads to the predicted Ising critical 
exponents. 

In the present paper, we want to consider the classical continuous spin model on 
introduced by Griffiths 0, an Ising model with spins taking values continuously between 
—1 and and prove that in this case magnetization and percolation transitions coincide. 
We also present a detailed numerical study of the model with simulations on the 2d square 
lattice. 

For SU(2) lattice gauge theory it was shown in the strong coupling limit [^j that the 
partition function can be written in a form which, apart from a factor which depends on 
the group measure, is the partition function of the classical spin model, with spins varying 
continuously in some bounded real interval 0]. These continuous spins thus are the natural 
counterparts of the Polyakov loops in SU(2) gauge theory. Following the pioneering work 0], 
the relationship between the thermodynamical features of spin models and the properties of 
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the corresponding geometrical clusters have in the past years received considerable attention 
pj]]). We want to address here the random percolation content of the continuous 



see 



spin Ising model |Q obtained through a Fortuin-Kasteleyn transformation of the partition 
function. A specific class of geometrical clusters shows a percolation transition whose 
critical behaviour matches the conventional Ising counterparts. In particular, it can be shown 
that up to finite constants, the spin magnetization equals the probability of long range cluster 
connectivity; moreover, the spin susceptibility becomes equal (up to the same constants) to 
some mean cluster size [|T^. 



The numerical method to be used here is the Wolff algorithm [0. It is based on the 



formulation of a single cluster update algorithm for spin systems with unprecedented per- 
formances regarding the problem of critical slowing down behaviour near phase transition 



12| . It turns out that the Wolff clusters exhibit critical percolation behaviour at the point 
of the Ising phase transition. Very recently, this specific property of Wolff clusters has been 
put into a rigorous framework [jlO|, leading to a concept of Wolff measures which provides 



the theoretical basis for the method to be used below. We shall extend this work to the 2-d 
continuous spin model and corroborate the results through numerical simulations. 

2 The Model. 

Let A C be a finite lattice with edge set E and vertex set V = {!,.., A^} and S a random 
configuration of spin variables S = {Si} C S^, where the Si^i^v take values independently 
in S = [—1, +1]. Consider the Hamiltonian 



HiS) = -J2 J^,ASj^ (2.1) 



<id> 



where the sum is over all nearest neighbour pairs (< i,j >) of spins and Jjj > 0. Rewriting 
the spin variables St = kai, where (jj = sign{Si) = ±1 and k = \Si\ E [0, 1], one has 

H{S) = H{a, l) = -J2 JiM^'i^'j (2-2) 

2.1 Joint Distribution over Spin Amplitudes and Bond Variables 

Following Edwards and Sokal [§], the random cluster representation (RCR) 0, ^ of this 
model is derived in the following way. Let n C {0, 1}^ be a configuration of bond variables 
on A where, 

- Uij = 1 means that there is a bond between sites i and j in V, 

- riij = is the opposite event. 

Consider the following joint probability distribution of the random variables {a, I, n), 



P(S,n) = P(Z,a,n) = Z^' n ^^"''^'''^ hA.,A.^, + (1 -p.,)5„,,,o] 

<i,j> 

= Z^' J] w^,j(Z,o-,n), (2.3) 



<i,j> 
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where pij = 1 — Q-~^p-'i,j'-ih is a weighting factor to be interpreted as the probabihty for a hnk 
Uij to exist {uij = 1) between two nearest neighbour spins Si and Sj that share the same 
sign {Sa^aj = 1), and (1 — Pij) the probabihty of having no hnk {riij = 0) between them. Za 
is the partition function of the model given by 

= J] / n ^^^(^*) J^M^^^J]^ (2-4) 

W} '' i <i,j> 

where dfi is the probabihty distribution of the spin amphtudes li and jS is the inverse tem- 
perature. 

Starting with the joint distribution ( |2.3| ), it is straightforward to express the weights 
Wi3j{l, cr, n) in the foUowing way 

wf},j{i,cT,n) = Yi e^-^'''^'^' n n -p^i'^- 

<i,j> <i,j> <i,j> 

nij=l nij=0 

Now, we write n for a bond configuration fulfiUing the compatibihty condition that a bond 
hij exists between sites i and j (hij = 1) iff ai = aj and caU c(n) the number of clusters 
of bonds in the configuration n. Then, the weights Wj3^j take the form (FK-representation) 

«^/3,j(z,n)= n e'^'-^^'M 2^^"^ n n (^-p^A- (2-5) 

<i,j> \ <i,j>€n <i,j><^n j 

These weights are called the Wolff weights of the random cluster distribution. The joint 
distribution on the configurations (Z, n) provide the right framework to relate the critical 
behavior in the original spin system and its associated percolation representation. 

Remark 2.1 The Wolff' weights Wi3j{l,fi) differ from the (FK) weights of the Ising case 
(see 1^^) by the term ^P-^i,i'-ih ifidi reflect the measure on the spin amplitudes in the 

<i,j> 

continuous case. 



3 Wolff Cluster Algorithm and Distributions 

Non-local cluster Monte Carlo (MC) algorithms have brought significant improvements in 
the simulation of Ising models near criticality. Starting from the ground-breaking work of 
Fortuin and Kasteleyn 0, which relates the partition function of spin systems with that of a 
correlated percolation model, Swendsen and Wang derived a non-local cluster algorithm 
which drastically reduces the critical slowing down phenomena near the transition point, 
with a dynamical critical exponent near 0.25, where z is defined by 

r = r. (3.1) 

Here r is the correlation time in MC simulations (measured in MC steps per site) and ^ is 
the correlation length. For a system of size L, near the critical point, ^ scales like in two 
dimensions. In the case of a local update algorithm like Metropolis or Heatbath, z is found 
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to be close to 2, so the required time to reach stable configurations at criticality in this case 
is of order L^, reducing noticeably the possible size of samples to study. 

Building on the Swendsen-Wang idea, Wolff |12| improved the method (see below) and 
derived a non-local update MC scheme with a dynamical critical exponent z smaller than or 
equal to Swendsen-Wang' s in any dimension and furthermore easier to implement. 



3.1 The Wolff Algorithm 

First let us recall briefly the main features of the Wolff method. Consider the 2-D ferromag- 
netic Ising Model on with coupling constant J and inverse temperature /3. Starting from 
a randomly chosen spin o"o, visit all nearest neighbours and 

- with probability pf, = 1 — exp{—2j3J), 

- and only if they have the same orientation as (Jq, (IH-l) 

include them in the same cluster as uo; spins not satisfying both conditions are excluded. 
Repeat iteratively this procedure with newly added spins in the cluster until no more neigh- 
bours fulflU the above compatibility condition (///.I). Now flip all the spins in that cluster 
with probability 1. After that, erase all the bonds and start this procedure again. It turns 
out that this dynamics verifles the detailed balance condition, i.e. it samples the Gibbs 
distribution of the Ising model (see The distinguishing feature of the Wolff method 

compared to Swendsen- Wang's is that, in the latter (following Fortuin and Kasteleyn 0) 
one needs to build, with the same growing probability pb as before, all possible clusters of 
like spins and then, with probability 1/2, flip all the spins in those clusters. Then all the 
bonds are erased and one starts again from the newly created spin conflguration. It can 
also be shown that this method verifles the detailed balance condition P|. In summary, the 
following remarks can be made: 

- The building procedure of an individual cluster in both methods is clearly the same, 
thus the Wolff cluster belongs to the set of Swendsen-Wang clusters. 

- When the Wolff cluster is built, the randomly chosen spin has obviously higher prob- 
ability to fall in a large Swendsen-Wang cluster than in a smaller one. 

It results that the distribution of Wolff clusters is given by the distribution of Swendsen- 
Wang clusters (see below) modifled by an additional weight that takes into account the size 
of the clusters. 

Remark 3.1 A (rather non economic) way of building the Wolff cluster would be to con- 
struct the set of all M Swendsen-Wang bond clusters Cf^ , i G {l,..,Af}, then to pick at 
random a lattice spin and flip with probability 1 all the spins belonging to the cluster 
C'^^((Tx). Obviously, the randomly chosen spin cTx has higher probability to sit in a large 
cluster. 

Remark 3.2 In the case of our model, the bond probability pb{Si, Sj) = 1 — exp{—2(3Jijlilj) 
takes also into account the spin amplitudes. It turns out that, the Wolff dynamics (as de- 
scribed above) is no longer ergodic (because spin amplitudes are locally conserved) and must 
be supplemented by a local update method (like Metropolis or Heatbath) in order to respect 
detailed balance. In practice it is sufficient to insert a Heatbath sweep after several Wolff 
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sweeps have been performed. This of course slightly alters the gain brought by the non-local 
cluster update method concerning critical slowing down near the phase transition. 

The proof that the Wolff algorithm applied to our infinite-spin model fulfills both ergod- 
icity criteria and detailed balance conditions will be omitted here since it follows closely the 



derivations that can be found in 15, 18, EQ, 21 



3.2 Conditional Distributions and Properties of Wolff Clusters 

We state in this section our main results concerning the properties of the RC distribution. 

First, it is straightforward to see that integrating P(Z,cr,n) over n gives the marginal 
distribution of the model with Hamiltonian (123 



Proposition 3.1 The marginal spin distribution in ( ^.Sj ) gives the Boltzmann weight asso- 
ciated with Hamiltonian \2.2^ ). 

Tracing over bond variables in (|2.3| ), one gets 

1 1 

nij=0 <jj>nij=0 

=zx' n {^''-■'''Hp.a.^. + - p^^)]} 

<i,j> 

using the expression for pij and re-arranging terms, this can be rewritten as 

<i,j> 
<i,j> 

thus we get for this marginal distribution 



• ^A ' 

<^,3> <«,J> <«,J> 

which is just the Boltzmann weight for the spin configuration S of our model. 

The nice features that make Wolff clusters useful reside in the following propositions: 

Proposition 3.2 The magnetization M of the spin system is equal to the probability of long 
range connectivity in the percolation model. 

Considering the joint bond-spin distribution (|2.3|) with weights (|2.5D, one can ask: given 



a bond configuration, what is the conditional distribution on the spins ? Let us assume 
boundary conditions for the spins (which amounts at fixing ZjCTj = 1 on the boundary). First, 
due to the compatibility condition (///.I), if site k is connected to the boundary through 
some path in the percolation model, then necessarily ak = +1. The same is true with (— ) 
boundary conditions due to the sign symmetry of the exponential factor in the statistical 
weight. Second, a site that is not connected to the boundary is equally likely to be in a 
±1 state. From this it follows that M is exactly given by the probability of connectivity to 
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the boundary, thence long range order in the spin system is the same as percolation in the 
random cluster model. 



A more detailed argument will be found in . It is actually easy to show that absence 
of percolation is equivalent to zero magnetization. Consider again the random cluster model 
(RCM) characterized by the joint distribution (|2.3| ). First let us state the 

Definition 3.1 Let PA,j(/5) be the probability for a site i to be connected to the boundary of 
a finite domain A C through some path of active bonds in the percolation model. We say 
that there is no percolation in the RCM iff 

V?:GA, lim PA,^(/?) = IPoo,^(/?) = 

A—flP 



Proposition 3.3 Consider the joint bond-spin distribution ( \2.^ and suppose that for the 
infinite-spin model at inverse temperature j3 there corresponds no percolation in the joint 
correlated geometrical bond representation; then the magnetization vanishes. 

The argument is as follows. Let cr„ = [ai^,ai^, ■■I'^in) ^ finite set of spins and Mj^\(T„) 
the corresponding magnetization in the spin model with boundary conditions (03). Let 
(53 — e), the boundary condition identical to (03) except that each spin value on the boundary 
has been lowered by an amount e. Then one has 

where cr„ + e = (crj^ + e, (Tjj + e, .., cTj^ + e). Clearly the geometrical clusters corresponding 
to these two boundary conditions are the same, thus for all bond configurations in which 
sites ^l,^2,..,^n are not connected to the boundary, the corresponding contributions to the 
magnetization are the same. 

Now, using the uniqueness of the limiting state at given it is easy to conclude that, 
in the infinite volume limit of definition ( |3.1[ ), absence of percolation is the same as zero 
magnetization. 

From the proposition above follow also interesting results concerning high temperature 
decay of correlations for functions of the spin variables when their relative distance goes to 
infinity but we will not enter this question here and refer to . 



The stronger statement of proposition ( p.2| ) whose details of proof will also be found in 
the last reference (see also [llO| , |Tl| ) is the following: not only the onset of percolation is the 



signature of the magnetic ordering but moreover, there exist a constant C depending only 
smoothly on /5 such that 

IPoo(/5) > M(/3) > CPoo(/3) (3.2) 

Explicitly, there is spontaneous magnetization if and only if there is percolation and the 
critical behaviour is the same for these two quantities. Namely, considering on one hand the 
order parameter exponent (3, defined for the spin system by M ~ (Tc — T)^, where Tc is the 
critical temperature, and on the other hand the order parameter for the percolation model, 
i.e. the fraction of sites F^o in the percolating cluster, which behaves like [p — Pc)^, where 
Pc is the percolation threshold. It follows from ( |3.2D that 13 = 13. 



A similar result as ( p.2|) can be proven (see |jT3|) between the magnetic susceptibility of 
the continuous spin system and its geometric counterpart. 

In the following we shall carry out lattice simulations to determine the critical expo- 
nents for cluster percolation, to illustrate that these indeed lead to the values of the Ising 
universality class. 
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4 Numerical results 

We have performed extensive simulations of our model, choosing six different lattice sizes, 
namely 64^, 96^, 128^, 160^, 240^ and SOO^. 

Our update step consisted of one heatbath for the spin amplitudes and three Wolff 
flippings for the signs, which turned out to be a good compromise to reduce sensibly the 
correlation of the data without making the move be too much time-consuming. Every five 
updates we measured the quantities of interest. The physical quantities are the energy and 
the magnetization of the system. The percolation quantities have been evaluated in this 
way: for a given configuration we consider the sign of the spins corresponding to the sign of 
the global magnetization, say up. Then we form clusters of spins up using the local bond 
weights 1 — exp{—2{J/kT)siSj). Once we embody all spins up in clusters we assign to each 
cluster a size s which is the number of spins belonging to it. We say that a cluster percolates 
if it spans the lattice in both directions, that is if it touches all four sides of the lattice. This 
choice was made to avoid the possibility that, due to the finite lattice size, one could find 
more than one percolating cluster, making ambiguous the evaluation of our variables. 

Finally we assign to the percolation probability P the value of the cluster size of the 
percolating cluster divided by the number of lattice sites (P = if there is no percolation) 
and to the average cluster size S the value of the expression 



5 = E 



where is the number of clusters of size s and the sums exclude the percolating cluster. 

After some preliminary scans of our program for several values of the temperature f3 
{(3 = J/kT), we focused on the beta range between 1.07 and 1.11, where the transition 
seems to take place. The number of iterations for each run goes from 20000 (for (3 values 
close to the extremes of the range) to 50000 (around the center of the range). We interpolated 
our data using the density of states method (DSM) p2| . 

To locate the critical point of the physical transition we used the Binder cumulants 

Figure 1 shows Qr as a function of (3 for the different lattice sizes we used. The lines cross 
remarkably at the same point, which suggests that also in our case Qr is a scaling function. 
As a numerical proof we replot our lines as a function of (/3 — PcntjL^^'^, choosing for the 
exponent u the Onsager value 1. The plot (Figure 2) shows that indeed Qr is a scaling 
function. 

To get the critical point of the percolation transition we used the method suggested 
in [^, which is based on the same principle. For each f5 value we count the number of 
configurations where we found a percolating cluster and we divide it by the total number 
of configurations. If we plot the results as a function of /3 for the different lattice sizes the 
corresponding lines should cross at the threshold of the geometrical transition (Fig. 3). 

As the figures show, the agreement between the two thresholds we found is excellent. 

For the evaluation of the exponents we used the method [|^. The results we got are 
reported in Table I. 
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Figure 1: Binder cumulants as a function of /3 for our six lattice sizes. 
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Figure 2: Rescaled Binder cumulants. We took (3crit = 1-0932 and for the exponent u the 
Onsager value 1. 
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Figure 3: Fraction of percolating samples as a function of /3 for our six lattice sizes. 
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The two thermal exponents' ratios are in good agreement with the Onsager values. Un- 
fortunately we encountered some troubles in deriving the percolation /S/v exponents' ratio 
because the best fit values change rapidly in the little range of beta values for which the 
corresponds to the 95% confidence level. Consequently, the relative error oi l3/v turns out to 
be large. For a more precise evaluation of this ratio much higher statistics or higher lattice 
sizes seem to be necessary. 

However, the value of the other exponents' ratio ^/v is in good agreement with the 
corresponding Onsager value 7/4 and its error seems to exclude the possibility for it to be 
instead the random percolation ratio 43/24. 



5 Conclusions 

For the continuous spin Ising model, the percolation transition of Fortuin-Kasteleyn clusters 
is indeed equivalent to the magnetic transition due to spontaneous Z2 symmetry breaking. 
We have also determined by means of lattice simulations the critical exponents of the perco- 
lation variables and shown that they belong to the Ising universality class. The simulation 
method used here can also be employed in the study of percolation in SU{N) gauge theo- 
ries, where analytic proofs do not exist. We expect that they will eventually allow a general 
demonstration that the deconfinement transition in such theories can indeed be defined as 
Polyakov cluster percolation. 
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